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ABSTRACT 

Utility (e.g., sum-rate) maximization for multiantenna broadcast and 
interference ciiannels (with one antenna at the receivers) is loiown to 
be in general a non-convex problem, if one limits the scope to linear 
(beamforming) strategies at transmitter and receivers. In this paper, 
it is shown that, under some standard assumptions, most notably that 
the utility function is decreasing with the interference levels at the 
receivers, a global optimal solution can be found with reduced com- 
plexity via a suitably designed branch-and-bound method. Although 
infeasible for real-time implementation, this procedure enables a 
non-heuristic and systematic assessment of suboptimal techniques. 
In addition to the global optimal scheme, a real-time suboptimal al- 
gorithm, which generalizes the well-known distributed pricing tech- 
niques, is also proposed. Finally, numerical results are provided that 
compare global optimal solutions with suboptimal (pricing) tech- 
niques for sum-rate maximization problems, affording insight into 
issues such as the robustness against bad initializations in real-time 
suboptimal strategies. 

Index Terms — Nonconvex optimization, branch-and-bound, 
interference channel, multiple-input single-output channel 

1. INTRODUCTION 

Precoding and power control are well studied strategies that sup- 
port high spectral efficiency in wireless network with multiple an- 
tenna transceivers, when channel state information (CSI) is available 
at the transmitters. Several system-wide objective functions have 
been considered in the literature for precoding and power control 
optimization of broadcast channels (BCs) and interference channels 
(ICs). Some of these problems are convex, for example power mini- 
mization 1 1 1 or SINR balancing for the multiple-input single-output 
(MISO) BC |2], and thus solvable with standard techniques in rea- 
sonable (polynomial) time. However, in general, the problems at 
hand are non-convex. Unlike convex problems, non-convex prob- 
lems typically do not afford efficient (i.e., polynomial-time) algo- 
rithms that are able to achieve global optimality |9|. For example, it 
is known that the weighted sum-rate maximization (WSRM) in par- 
allel IC channels, where interference from other users is treated as 
noise (a non-convex problem) is NP-hard [3 1 (this result extends also 
to BC as a special case). 

In this paper, we address the global minimization of a system- 
wide, in general non-convex, cost function with respect to the trans- 
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mit covariance matrices {Qfc}. Among global techniques, branch- 
and-bound (BB) algorithms are methods to solve general non-convex 
problems |5 |, producing an e-suboptimal feasible point. BB meth- 
ods have been already introduced to solve non-convex power control 
problems, although so far only multi-user single-input single-output 
systems have been addressed in f6l and references therein. 

In this paper, we propose a novel BB framework for global opti- 
mization of a problem formulation that includes, for instance, MISO 
BC and IC WSRM with general convex power constraint. The pro- 
posed BB approach is based on the observation that a fairly general 
set of cost functions that arise in communication's problems, albeit 
non-convex, possess a Partly Convex-Monotone |7| structure. This 
structure is satisfied whenever one can identify a suitable set of inter- 
ference functions f, ({Qfc}), for which the following hold: (;) The 
cost function is convex in the transmit covariance matrices {Qfc} 
once the interference functions fi ({Qfe}) are fixed; (n) The cost 
function is monotone in the interference functions fi ({Qfe}). We 
design the BB scheme to exploit the Partly Convex-Monotone struc- 
ture of the problem. Branching is performed in a reduced space (of 
the size of the set of all feasible interference level vectors fi ({Qfe})), 
instead of the original feasible space (of the size of the set of all fea- 
sible covariance matrices {Qfe}). Bounding is efficiently carried out 
by solving only convex optimization problems. 

In addition to the reduced-space BB method, we propose a sub- 
optimal algorithm that attains quasi-optimal performance with poly- 
nomial complexity. This algorithm reduces to the distributed pricing 
scheme of (4|, when applied to sum-rate maximization problems. 
Numerical results are provided to compare the global optimal solu- 
tion based on BB, the suboptimal (pricing) technique and the non- 
linear dirty-paper coding scheme. 

Notation: The Boldface is used to denote matrices (uppercase) 
and vectors (lowercase); (•)"^ and (•)^ denote the transpose and the 
Hermitian transpose, respectively; Tr (■) denotes the trace of a ma- 
trix; E [■] denotes the expectation operator. Moreover, given a vec- 
tor X we address its Z-th component as [x]p and the vector inequality 
X =^ y means that [x]j < [y]; Vi. Finally, unless otherwise specified, 
we address the set of covariance matrices {Qfe}^^ as {Qfe}. 



2. SYSTEM MODEL 

We model a multi-user communication system consisting of K 
transmitter-receiver pairs (or users). The fc-th user has Nk transmit 
antennas and one receive antenna (MISO system). The signal at the 



k-th receiver is given by 



2.2. Examples 



signal 



(1) 



inteiference 



where G c^*"^^ is tiie fc-th transmitter's signal, hfcj G C^^^*" 
accounts for the channel response of the MISO link between the 
fc-th transmitter and the j-th receiver, and Wk G C'^'^^ mod- 
els the additive white Gaussian noise (AWGN) at fc-th receiver: 
Wk ~ CN'{Q,(j\). Assuming capacity-achieving Gaussian code- 
books, we define the coiTelation matrix of the fc-th transmitted signal 
as Qfc = E [xfeX^] . While model lUll accounts for an IC, a BC can 
be obtained as a special case by setting \\jk = Vj, k. 



2.1. Problem Formulation 

Due to multi-user interference, the system performance depends on 
the transmission strategy of every user, i.e., on the set of covariance 
matrices {Qfc}. We consider the minimization with respect to {Qfe} 
of a system-wide cost function / (to be defined below) under a gen- 
eral convex set constraints Q: 



miri /({Qfe},fi({Qfe})) 



(2) 



By defining a set of L auxiliary variables i, problem ^ can be recast 
in the equivalent form 



(P) min /({Qa,i) 

{Qfc}6S,i 

s.t. i = fi ({Qfc}) 



(3) 



The equivalence means that if {Q^} is a solution to l|2j, then 
({Qfc} ,fi ({Qfc})) is a solution to ([3}. Conversely, if ({Q^} 
is a solution to l[3]l, then {Qfc} is a solution to l|2j. 
We further make the following assumptions: 

Al The L interference levels are given by the real vector function 
fi ({Qfc})> ojfine with respect to {Qfc}, that is bounded in 
the L-dimensional rectangle [i^in jmaxj ^ t^l ^j^^ 

th component satisfies if''" < [ft ({Qfc})], < ^1""" for I = 
!,...,!/). For instance, we typically have L = K and the 
interference level at the fc-th receiver reads [ft ({Qfc})]^. = 

Tlf=i,j^k hjfcQjhj^; 

A2 The cost function / ({Qfc } , i) is a real scalar function that is: 
continuous in ({Qfc} , i); monotonic increfljin^ with respect 
to i G [i™'", i ] for fixed {Qfc} G Q; convex with respect 
to {Qfc} for fixed i G [i""'", i™""] ; 

A3 The set Q is closed and convex. For example, Q may 
be the set of positive semidefinite covariance matrices 
{Qfc )p 0} satisfying the generalized power constraints 
Y^k^^ Tr (Afc,fQfc) < P, for ^ = 1, . . . , D, where {Afc,,} 
are positive semidefinite matrices (possibly A.k^t = if 
fc-th user doesn't belong to ^-th constraint) and {Pi} are 
non-negative coefficients. This definition includes some 
important special cases studied in the literature, such as per- 
antenna, per-group of antennas, the classical sum-power or 
the interference constraints in cognitive radio scenarios. 

Throughout the paper, we refer to problem (|3] as (P). We next 
provide examples of problems that satisfy these assumptions. 



An example of cost function included in our framework is the a- 
faimess criterion 18]: /({Qfc}) = -Wfc/a ('"fc ({Qfc})), 

where Wk is a positive constant, fa is an increasing strictly concave 
function defined as 



U (r) : = 



log r 
(1-a) 



if a = 1 
otherwise 



(4) 



and rfc ({Qfc}) is the fc-th user's rate, which depends on covariance 
matrices {Qfc} and on the channel scenario. The a-faimess criterion 
reduces, as special cases, to the WSRM problem (a = 0) or the 
proportional fairness problem (a = 1). Moreover, as a becomes 
large, it converges to the max-min fairness problem |8]. 
In the following we present some examples of channel scenario that 
can be addressed within our framework: 

• Parallel MISO IC: The fc-th transmitter operates over Lc par- 
allel subcarriers, it has power constraint Pk and has knowl- 
edge of channels {\\jki} for j = 1, . . . , K and VZ. 
The minimization of the (p, a) -fairness cost function reads 



min -Wkfa (rk) 



s.t. 







rfc<Elog 1+- „K 

i = l V ''fcl+Z^j^l.jT^fchjfclQjihffc, 

Tr (Efii Qki) < Pk 



(5) 

Vfc 
Vfc 



Defining [ft ({Qfci})],+£^{fc_i) = E^i j^^fc hjfciQjihj^, 
Vfc, UndQ= |Qfc, ^ Vfc, Z I Tr (j^ffi Qki) < Pk Vfcj, 
problem (|5} is recast into (P). Also, ft ({Qfc;}) G [O.i'"'"'] 
where i™''" is a proper upper bound on interference, always 
available since Q is bounded (finite power constraints). 

Parallel MISO BC: This scenario is obtained from (|5} by set- 
ting hjki ~ hfcj Vj and imposing a sum-power constraint 



By suitably modifying the same arguments, the proposed framework can J non-overlapping sub-rectangles (i.e., -^t^' 



3. PROBLEM SOLUTION VIA BRANCH-AND-BOUND 

In this section we show that, adopting standard BB techniques (see 
(5] and |7|), problem (P) can be optimally solved by means of 
an efficient BB that exploits the structure dictated by assumptions 
(A1-A3). The BB algorithm is fully characterized by two proce- 
dures: branching and bounding. These are iteratively performed un- 
til the solution's suboptimality falls below some prescribed accuracy 
e. In the following we explicitly tailor those procedures to prob- 
lem (P) and we show convergence of the proposed BB algorithm to 
the global optimal solution of (P). For readability's sake, we define 
Q : = {Qfc jEl and we address an interval as : — [a, b] , meaning 
that c G 7W ^ [a], < [c], < [b], for Z = 1, . . . , L. 

3.1. Branching Procedure 

A partition set Vt of rectangles {M} in the space R^, each labeled 
with a lower £(8 i-M) and upper il® (M) bounds, is given. By 
splitting a rectangle that satisfies Mt G argmiuAigPt £» i-M) in 

= 0and 



handle an analogous but more general case where / ({Qfc} 
monotone increasing in 1+ and monotone decreasing in i~ 



results 



^Here and in the following, the expression Q G Q stands for {Qfc} G Q. 



U/=i -^F = -Mt), the enhanced partition Vt+i = {VtXMt} U 
I^Wtl is obtained. Lower and upper bounds for each sub-rectangle 
in |7V!f I are then obtained via the following bounding procedure. 

3.2. Bounding Procedure 

Exploiting the Partly Convex-Monotone structure of problem (P), 
for every rectangle M = j^jmm jmax-j ^ Vt, & lower bound 
£b {M) is evaluated by solving the following problem: 



4. SUBOPTIMAL SOLUTION 

While the proposed BB algorithm always converges to the global op- 
timal solution and has reduced complexity with respect to a general- 
purpose implementation of BB, it is still feasible only for offline 
simulation. In this section we propose a suboptimal algorithm with 
polynomial complexity that extends the distributed pricing schemes 
of | 4 | to the more general class of problem (P). 
Exploiting the Partly Convex-Monotone structure, problem (|3) can 
be equivalently reformulated as the non-convex problem: 



: = min/ Q,r" 
QeQ V 

I •mill , n //-v\ , 

s.t. 1 =^ fi (Q) =^ 1 



(6) 



Thanks to assumptions (A1-A3) two fundamental results can be ver- 
ified: (0 problem ([6j is convex since the cost function / (Q, i) is 
convex for a fixed i and the constraints form a convex set, {ii) using 
standard convex optimization arguments, it can be shown that this 
bounding procedure satisfies the natural condition: 

M' (Z MUSI'S {M') > £<8 {M) . (7) 

Moreover, denoting with Q'~'^ ' the optimal solution of problem (|6), 
a valid upper bound Haj (A^) is obtained by evaluating the function 

atQ'~»),i.e.,il<B (M) : = / (^Q'~'='',fi (q*~'^')). 
Finally, the algorithm checks if the prescribed accuracy is met (i.e., 
if miriil<B ~ min£<B (M) < e) otherwise it goes back to the 
branching procedure. 

3.3. Convergence Analysis 

Here we proves convergence of the proposed BB algorithm. 

Lemma 1. The proposed BB algorithm (which is performed in the 
reduced space spanned by interference levels/variable i }, is conver- 
gent to a global optimal solution of problem (P). 

Proof. As explained above, since the chosen bounding procedure 
satisfyies O, the BB algorithm generates a sequences of partition 
sets {A^t} collapsing to a point Ht^oo ~ ^* (recall that Mt 
is the rectangle selected for splitting at the i-th branching iteration). 
In order to prove convergence we need to show that, as the size of 
rectangle Mt gets smaller, 11(b (Alt) — £58 [Mt] is also sufficiently 
small. The proof follows standard arguments |5 |. This is shown in 
Appendix. □ 

3.4. Broadcast WSRM Example 

Considering the BC WSRM scenario (see Sec l2.2t . for a given inter- 
val M 



mm jmax] 



,1 



the evaluation of a lower bound £58 results 



in the following convex problem 



£« {M) 



mm 



'Wk log 1 + 



(8) 



s.t. 



Defining {Qfe} the optimal solution of {HJ, a valid upper bound is 
given by JI95 (M) : = EfLi -Wfc log where it = 



mm sup 



mm 
qeQ 



[/ (Q, 



A^fi (Q) 



A 1 



(9) 



where A is the Lagrange multiplier associated to the affine constraint 
i = fi (Q). Building on ([9}, in the following table we formalize the 
proposed suboptimal algorithm. 

Algorithm 1 - Suboptimally solve problem (P) 

0: Set e\,e\ 

1: Initialize A = A 

2: Initialize i — i 

3: Evaluate Q* — arg rnin/ ^Q, i^ + A fi (Q) 

4: If ||T-fi(Q*)|| >ei 

5: Updatel ^ fi [Q") 

6: Go back to step 3 

7: elself A - ^^A^l ^ > £x 



Update A : 



9: Go back to step 2 
10: end 



Q=Q*,i=i 
9/(Q,i) I 



Since a stationary point of this algorithm fulfills the necessary 
Karush-Kuhn-Tucker (KKT) conditions of problem {3]), if the algo- 
rithm converges, it attains a local optimal point of problem ([3]l. 
It is worth noticing that, by specializing our framework to the case 
when the cost function / (Q, i) is the WSRM (i.e., fa ifk) ~ 
Wk in {5}) the Lagrangian multiplier A plays the role of the inter- 
ference prices defined in the distributed pricing algorithm f4\. Thus 
Algorithm [T] can be seen as a generalization of distributed pricing 
technique with an arbitrarily cost function and arbitrary interference 
functions (satisfying assumptions A1-A3). 

Finally, since the problem at hand is non-convex, initialization of 
the parameters A and i results crucial for performances and con- 
vergence. In Sec|5]we assess the performances of this technique in 
relation to the global optimal solution evaluated via BB algorithm. 



5. NUMERICAL RESULTS 

We assess the perfomance of the two proposed techniques: a multi- 
user linear precoder optimized via (;) the efficient BB algorithm (BB 
- LB and UB); (ii) the suboptimal Algorithm[T] 
We consider the sum-rate utility function (i.e., in (|5), fa (rk) = fk 
Vfc and Wk = 1 Vfc). In BB algorithm, the solution's accuracy is 
e = 10""^, while, in Algorithm[T] we run two different price initial- 
izations (Afe = 10~^ Vfc and Xk — I Vfc) and for both we initialize 
ik = u'i Vfc, selecting (j| = 1 at each receive antenna. 
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Algorithm 1 Algorithm 1 

(lnit:\=1Vk) (lnlt:\= 10"^ Vk) 
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Fig. 1. Sum-Rate versus transmitting power for the single-carrier 
(Lc ~ 1) BC scenario. The figure compares the two proposed lin- 
ear precoding algorithms: the optimal BB algorithm (BB - LB and 
UB) and the suboptimal Algorithm [T] considering two different ini- 
tializations (Afc — 10^^ Vfc and Afc = 1 Vfc). The optimal non-linear 
DPC technique is also plotted as a reference. 



Fig[T]shows the sum-rate versus the transmitting power for a single- 
earner (Lc = 1) BC channeO where a A'^ = 4 transmit anten- 
nas base-station serves K = A single-antennas users, subject to a 

sum-power constraint, Tr ^X^fLi Qfe) — Ptot- The sum-capacity 
achieving non-linear technique Dirty Paper Coding (DPC) is also 
plotted as a reference. 

It can be noticed that the suboptimal Algorithm [T] while showing 
near-optimal performance at several power levels, happens to be 
quite sensitive to initialization. For instance, initialization Afc = 
10~^ Vfc yields a suboptimal slope in high power regime, as ob- 

= 2MB both initial- 



served for Ptot 



> 31, and, at Ptot 



izations lead to highly suboptimal performances. A last observation 
pertains to the significant gains of non-linear DPC with respect to 
linear precoding at high power regime. 

Finally, not to confuse the reader, since a utility function (sum-rate) 
instead of a cost function is plotted, in fig[Il the lower bound results 
as the maximum feasible value while the upper bound is the maxi- 
mum upper bound among BB partitions. 



6. CONCLUSIONS 

This work presents a global optimization framework for the mini- 
mization of non-convex cost functions in MISO BC and IC channels. 
Examples are given for the general a-faimess optimization consider- 
ing parallel IC and BC channels. Knowing the global optimal solu- 
tion, even if impractical for real-time implementation, allows to as- 
sess the quality and to fine-tune (e.g., initialize) suboptimal schemes. 
In addition to the global optimal BB, we have proposed a real-time, 
hence suboptimal, algorithm that generalizes the pricing scheme of 
|4| . Extensions to MIMO networks are the subject of future work. 



^Due to space limitation, Fig[T] channels realization is available at: 

http : / / web . n j it . edu/ ~mr227 / papers / paper _BB_H_BC . mat 



7. APPENDIX 

We need to prove that, as the maximum length of the edges of A^t, 
denoted by size{A4t), goes to zero, the difference between upper 
and lower bounds uniformly converges to zero, i.e., Ve > 35 > 
OVMt C Mo size{Mt) < 5 =^ U'siMt) - S-'s{Mt) < e. 
For eachT £ Mt = [i™'", i™'"'] , we define the function 5 as 
the result of the following constraint optimization problem: 



mm 
QeQ 



s.t. i" 



/ (Q, i) 
=^ fi (Q) =^ i" 



Using this notation, the lower bound in (|6) is given by £<8 = 
5 (i""") , while an upper bound is given by its = 5 (i™'"') . 

From jointly-continuity of the function / (^Q)i^ with respect to 
^Q, (assumption A2) and from the definition of 5 , we have 
that 5 ^i^ is continuous in the norm of i (i.e., ||^j|)- It follows that 
also £<B and U<8 will result continue in j |Tj | , thus it holds 



Ve3(5 



concluding the proof. 
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